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The analysis of two-particle azimuthal correlation function in p + p and Au + Au 
collisions at y / ijvjv= 200 GeV is presented. The mean jet fragmentation transverse 
, momentum ^/{jjj and mean parton transverse momentum {k 2 ^) are inferred 

f^*) ■ from p+p collisions. The analysis of Au+Au data, where elliptic flow is a dominant 

co : source of two-particles correlation is also presented. We derive the upper limit on 

partonic azimuthal anisotropy based on jet-quenching in a Glauber model. We 
^ , show that some new mechanism, generating the azimuthal anisotropy observed in 

1} ■ final state, is needed. 
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1 Introduction 



With the beginning of RHIC operation heavy-ion physics entered a new regime, where 
■ pQCD phenomena can be fully explored. High-energy partons materializing into hadronic 

jets can be used as sensitive probes of the early stage of heavy ion collisions. Measure- 
ments carried out during the first two years of RHIC operation at v / s]v7V=130 and 
200 GeV exhibit many surprising features. The high-pj_ particle yield was found to 
be strongly suppressed in Au + Au central collisions [1, 2], and the particle distribu- 
tion in azimuthal space reveals large asymmetries attributed to sizable elliptic flow [3]. 
Furthermore, the azimuthal anisotropy in the region of p±>2.5-3 GeV/c, where the 
contribution from hydrodynamics is expected to be negligible, is found to be signif- 
icant and p±_ independent [4]. Another unique feature, found in RHIC data, is the 
observation of back-to-back jet correlation disappearance in central Au + Au collisions 
[5] . Although tremendous progress in the understanding of these phenomena has been 
made, an unambiguous picture merging all observables in a consistent way still has to 
be developed. 

The QCD medium effect on high-px particle production in heavy ion collisions 
is quantified with the nuclear modification factor Raa(p±) given by the ratio of the 
measured AA invariant yields to the NN collision scaled p + p invariant yields (eq. 
(1) in [2]). With no nuclear modification of particle yield, Raa(p±) = 1- The large 
decrease of Raa(p±) observed in central Au + Au collisions is commonly attributed 
to the multiple interactions of highly energetic quarks with excited nuclear medium. 
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Induced soft-gluon radiation carries off a significant fraction of each parton's energy, 
which leads to the depletion of the high-pj_ particle yield [6] . 

It is expected that the same mechanism of parton energy loss in the nuclear medium 
should also generate azimuthal anisotropies observed in the high-pj^ region. However, as 
it was first pointed out by E.Shuryak [7] and B. Muller [8], the influence of the nuclear 
medium on parton trajectory, although sufficient for excessive yield suppression, is not 
capable to generate azimuthal anisotropies larger than 10%. 

It has been suggested that the partonic "collective" hadronization might amplify 
the hadronic flow to the values of V2, as seen in experiment, if one assumes some sort 
of partonic coalescence [9, 10, 11]. The v-i values of identified 7r ± , K ± and pp [12], and 
Kg, AA [13] scaled by its quark content are apparently suggesting some kind of partonic 
collectivity. However, more data, especially identified particles at high-pj_, arc needed. 

Relatively good quantitative description of the high-p^ V2(p±) data was achieved in 
the model [14, 15] based on extension of the original "Color Glass Condensate" (CGC) 
model [16, 17]. The mono-jet (one gluon) radiation in classical CGC has been extended 
by including two-gluon production amplitude. This two-gluon radiation term generates 
two-particle azimuthal correlations, which are evidently strong enough to account for 
the measured strength of V2(p±)- The elliptic flow in this model does not respect the 
reaction plane, and the azimuthal anisotropy is solely of the two-particle correlations 
nature. However, the absence of the correlation between spatial distribution of particle 
emission and the orientation of the reaction plane in this model seems to contradict 
some of the other observation, e.g. reaction plane dependent HBT [18]. Moreover, the 
recent results from d+Au analysis [19] leave only small room for initial state phenomena 
required by CGC model. 

In the next section we will present the analysis of the two-particle correlation func- 
tion in p + p collisions, where the elliptic flow correlation is not present and the jet 
properties can be studied in a "clean" environment. 

2 Jet fragmentation in p + p collisions 

at v /s^=200 GeV 

Jets are produced in the hard scattering of two partons. The schematic view of such an 
event in the plane perpendicular to the beam axis is shown on Fig. 1. Two scattered 
partons propagate back-to-back from the collision point and fragment into the jet-like 
spray of final state particles (only one fragment of each parton is shown) . These particles 
have a transverse momentum Jj_ with respect to the partonic transverse momentum. 
The magnitude of jj_ has been observed to be p± independent (see [20] ) . 

In the case of collincar partonic collisions the two emerging partons would have the 
same magnitude of transverse momenta pointing to the opposite directions. However, 
the partons are carrying initial intrinsic momenta k± before the collision happens. This 
affects the outgoing direction and magnitude of partonic p± . It results in a momentum 
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imbalance (the partons' p± are not equal) and an acoplanarity (transverse momentum 
of one jet does not lie in the plane determined by the transverse momentum of the 
second jet and the beam axes). The jets arc not collinear and have a net transverse 
momentum yf (pj_) air = v2 • \J (kj_) of magnitude [22]. If the jet axes were known, 



parton 




, parton 



Pout 



Fig. 1: The schematic view of the hard scattering event. Two back-to-back scattered par- 
tons fragment producing high-px particles labeled as p xTrigger and p±, j± denotes the jet 
fragmentation transverse momentum and p± pa i r the parton pair's net transverse momentum. 

then the average projection of jj_ into the plane perpendicular to the beam axes can be 
determined in the following way 

(Lu*l) = (pj_)sin(|0i - <f> jet \) (1) 

where (p±) is the average p± of the particles produced in the fragmentation process. 
Here (j)jet and <j>i are azimuthal angles of the parton and i-th fragment respectively. 
However, we can only measure the relative angular dispersion between two jet fragments 
i and j which is <r N = y/((<fii - 4>j) 2 ) = a/2 y/( (^ - (f> je t) 2 ) = ■ (\<f>i - <t>jet\}- The 
final formula for (|j'j_j,|) is 

(\j ±y \) = {p x )-sin^ (2) 
V 71 " 

In order to extract y/ (k'j_) from measured width of away side correlation peak, ap, 
we used the equation for the average transverse momentum (IfWtl) out of the plane 
defined by prrigger and the beam axis (see [20] ) . 

(\Pout\) 2 = (\j±y\) 2 + x E ({\3±y\) 2 + mk±y\) 2 ) (3) 

where 

P± ■ P±Triqqer , 
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with p± being the momentum of the particles produced in the opposite direction to the 
P-LTrigger- For the case where the magnitudes of p ^Trigger and p± are similar, then 



(5) 



x E w -cos((|A0|)) = - cos( \ -a F ) 



The magnitude of (\p ut\) is expressed as 



(bout |) = (pj_)sin((|A0|)) = (px>sin(\/-£7i?) 
From (3) and using (2) one can derive the relation 



(\k± y \) = (P±) COS^z 



- tan z 

2 \ V 7T 



2 ' 



tan z | \/ — (Tjv 

7T 



(6) 



(7) 



In the limit of small angles, erjv — > 0, o> — > and equations (2) and (7) simplify to 

(p±) ■ a N 



(P-l) ' \Af - CT 



(8) 
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Fig. 2: Two charged hadrons area-normalized correlation function extracted from p + p data 
for both particles in 1.0< p± <1.5 GeV/c (left panel) and 2.5< p± <5.0 GeV/c (right panel). 
The solid line correspond to the result of the fit (9). 
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An example of the two-particle correlation function (CF), denned as a ratio of "real" 
and "mixed" particle pairs distributions C(A<p) — d ^£ l / dN ,[&^ d ; inp+p collisions is 
shown in Fig. 2 for (p±) = 1.2 GeV/c (left panel) and (p±) — 3.1 GeV/c (right panel). 
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Fig. 3: Measured width of the near angle peak crjv (upper left panel) and far angle peak a F 
(lower left panel) as a function of p±. Extracted values of (|i± y |) at given px-bins are shown 
on upper right panel. The dashed line corresponds to the truncated "Seagull" fit (see text). 
The lower right panel shows the (|fc_Li/|) with zeroth order polynomial fit (dashed line). 



The correlation functions are area-normalized (/" C(A(j)) = 2n). A significant ex- 
cess of small angle (A(p—Q, intra-jet correlations) and large angle (A<fi — ±ir, inter-jet 
correlations) is evident. The data were fitted by 

Fit(C,Y N ,<j N ,Y F ,<j F ) = C + Y N ■ Gauss{A(j) = 0,a N ) + Y F -Gauss(A<j) = ±n,cr F ) (9) 
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function, where C is the constant, Yjv, ajy, Yp and <7f are the yields and peak widths 
of near and far angle peaks respectively. 

The analysis of the correlation functions for various pj^-bins were performed. Ex- 
tracted values of o n and of are shown on Fig. 3. The angular width of the near and far 
angle peak is decreasing with pj_ as expected in case of jet fragmentation. The values 
of and were extracted. In case of (upper right panel of Fig. 3) 

one can observe the reduction of extracted for lowest pj_-bins. This is caused by 

p± truncation, sometimes referred as "Seagull effect" [23] (the j± at given p± cannot 
be larger then p± itself, and thus the measured ct/v and j± at p± ~ j± is accordingly 
reduced.) One can easily account for such truncation and the result of the fit procedure 
is shown as a dashed line on upper right panel of Fig. 3. 
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Fig. 4: The comparison of PHENIX preliminary measurement of (|jxs|) at \/sjvjv= 200 GeV 
with CCOR measurement done at y / ijvjv= 31,45 and 62 GeV (left panel). The comparison of 
V^pS pair ex tracted from PHENIX preliminary measurement on (\k± v \) according (10) shown 
on right panel. 



The average value is 

(|i ±3/ |) - 373 ± 16MeV/c. 

The average was extracted by fitting a constant to the data above 1.5 GeV/c. 

The average value found is 

(\k± y \) = 725 ± 34 MeV/c. 
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The comparison of with results of CCOR collaboration [20] (left panel on Fig. 4) 

shows a good agreement. The (|fcj_ y |) parameter could be compared to average RMS of 
the pair's momentum (pi) pair i which has been studied onp+p Tcvatron collider at 
Fermilab [24]. The pair's p± is related to the individual parton ( | fcj_ y | ) by the expression 

J^) . =V2-J<rt) = V^-{\k ±y \) (10) 

V pair V 

Extracted value of \J (pj_) ir =1.82±0.85 GeV/c is in good agreement with p+p data. 
This value could be also compared to (pj_)=1.80 ± 0.23 (stat) ± 0.16 (sys) GeV/c of 
PHENIX J/* measurement [21]. 



3 Two-particle correlation in AA 

Since the colliding heavy nuclei have a finite radius, the collision zone, except in very 
central collisions, has typical "almond" like shape with transverse dimensions of order of 
few fm. This geometrical azimuthal asymmetry is usually characterized by eccentricity 

<» - <y lrl > (") 

w < y 2 + x 2 > y ' 

where x and y are the coordinates of colliding nuclear constituents in transverse plane 
and b is the magnitude of the impact parameter. In the hydrodynamical approach [25], 
the original geometrical asymmetry is transformed into momentum space by means 
of multiple interactions between particles. The anisotropy in the momentum space is 
parametrized in the same way as spatial anisotropy 

< v 2 — v 2 > 

where v-i (p± , b) corresponds to the second harmonic coefficient of Fourier expansion of 
the particles azimuthal distribution 

^-f a [1 + 2v 2 ■ cos(2(0 - *))] (13) 

Here (f> is the relative azimuthal angle of the particle emission with respect to the 
reaction plane "J. In the case of "pure" flow, when the particles' correlation is completely 
generated by (13), the correlation function has a form 

dN f r . , dN f dN 



dAcj) 



C dNf dNf , o / a \ \ 



It is obvious that the azimuthal anisotropy v 2 generated by nuclear geometry cannot 
exceed the magnitude of eccentricity e at given impact parameter. 
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In the "sharp sphere" approach, where the nuclear density is assumed to be constant 
inside and zero outside sphere of radius R, the eccentricity of a collision at impact 
parameter b = \b\ is e = b/2R, and the geometrical vi could be as large as 1; but 
this approach is too idealized. When the more realistic nuclear density distribution of 
"Saxon- Woods" type is taken into account [26] , the maximum value of e barely reaches 
50% (see Fig.5). 




Fig. 5: Left panel: the N co u and N par t distribution with impact parameter b computed for 
Au + Au collisions at ^/sjvjv=200 GeV. Total inelastic cross section <7t ot =40mb, Au radius 
J?=6.38fm and nuclear density of 0.17 GeV/fm 3 were assumed. Right panel: The eccentricity 
dependence on b for N co u and N par t distributions calculated with Saxon- Woods and sharp 
sphere nuclear density distributions. The geometrical eccentricity b/2R shown with dashed- 
dotted line. 



However, the momentum anisotropy generated by collision geometry, irrespective of 
details of collision dynamics, can be only smaller. Since the anisotropy in the high-pj_ 
region cannot be generated by hydrodynamic flow (we will ignore the initial state effects 
described in [14]), the expected source of observed anisotropy is the process of energy 
loss of high energy partons in excited nuclear medium. The upper limit on w 2 generated 
by jet-quenching mechanism can be computed in a similar way described in [7]. 

The parton born in hard scattering events inside the collision almond has to travel 
the distance L in the medium before it reaches the outer space. The "primordial" 
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azimuthal asymmetry of the medium, where the hard scattering event happened, grad- 
ually vanishes as the medium expands. Since we are searching for an upper limit on V2 
generated by jet-quenching mechanism we can neglect any dynamic and expansion of 
the collision zone. 

The energy loss of the parton propagating through the medium could be charac- 
terized by absorption coefficient k (dimension fm -3 ). The probability of escaping the 
medium for parton born at the space point xq , which propagates along the direction n, 
could be expressed as 

/•OC 

f(x*o, n) = exp[— k / ds L-(xq + sri) ■ L + (xq + sn)]. (15) 
Jo 

(see Fig. 7). Here L±(x, y) = 2[R 2 — y 2 — (x ± 6/2) 2 ] 1 / 2 , nuclear thicknesses along the 
beam axis in case of sharp sphere approach, characterizes the nuclear density at space 
point (x,y) and it is proportional to the probability of hard scattering [7]. 




Fig. 6: Left panel: the relative absorption rates for various values of parameter k compared to 
the PHENIX tt° Raa for p ± >4 GeV/c [2]. Right panel: calculated «2 as a function of N par t 
for the same values of parameter n compared with the PHENIX preliminary data [3]. The 
solid line represents the limit from sharp sphere with k — > oo (see (17)). 



In more realistic Glauber model one has to replace L_L + by N part (x) distribution, 
which characterizes the density of excited nuclear medium. The probability distribu- 
tion of hard scattering is approximated with the product of two thickness functions 
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Ta(x, y)Ts(x, y) (see [26]). As it was mentioned, this approach provides only an up- 
per limit of V2 from the jet quenching mechanism. Any more realistic scenario with a 
detailed calculation of partonic interactions with expanding QCD medium could lead 
only to reduction of the resulting v 2 - 

The k parameter is the only unknown parameter in this calculation and could be 
estimated from the comparison of the absorption rate dependency on collision centrality 
expressed in terms of N part (see Fig. 6). Adjusting the absorption coefficient k compar- 
ing to PHENIX 7T° Raa for p± >4 GeV/c, we found the k=0.1 fm -3 provides the best 
agreement with the data. But the maximum v 2 in this case is 6%. For less realistic 
values of k the maximum value of v 2 reaches 10%. However, the values of v 2 (p±) mea- 
sured by the PHENIX collaboration using two-particle correlation techniques exceeds 
this maximum "jet-quenching" v 2 (p±) by a factor of 3. 




Fig. 7: The schematic view of particle production in nuclear collisions. Particle produced 
inside the overlap region propagates in the direction ft. 

The solid line on the right panel of Fig. 6 represents an idealized limit in sharp 
sphere approach, when k — > oo, and only radiation emitted from the surface is not 
absorbed [27]. The escape probability P osc is in the a — > oo limit 

P osc = / dhexp(— pah/cos9) oc cos^lo-^oo. (16) 
Jo 

where p stands for nuclear density, a for absorption cross section and the product pa 
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has similar meaning as k parameter. The v 2 can be obtained from integration 



/rOim pTt/2 
d(j)P esc ■ cos(2</>) = da d6>cos 



cos 1(9 + a) sin2a n 



tt/2 2 ' 2a 6a m 

(17) 

where a m = ^/R 2 - (b/2) 2 . 

As one can see on Fig. 6, even the unphysical model of black body radiation providing 
the upper "geometrical" limit on v 2 , still does not generate enough azimuthal anisotropy 
as observed in data. It is evident, that in order to to explain measured azimuthal 
anisotropics, some new underlying mechanism, other than just the "jet-quenching" , has 
to be invoked. 



4 Summary 

We have presented the detailed shape analysis of the two-hadron correlation function 
in p + p collisions at GeV. The average jet-transverse momentum (|j_Lj,|} = 

373 ± 16MeV/c was measured. The average partons transverse momentum was 
extracted by fitting the constant to the data above 1.5 GeV/c, and the average value 
of = 725 ± 34 MeV/c was found. The comparison of (|j_L y |) with measurement 

done by CCOR experiment [20] shows good agreement and no significant trend with 
increasing ^/JWn- The value of was compared to the values extracted from data 

take at Tevatron collider, and no significant deviation from the overall trend was found. 
The indirect measure of partonic primordial transverse momentum is the average p± 
of J/^> mesons. Such a measurement has been done by PHENIX collaboration as well, 
and the value of (p_i_)=1.80 ± 0.23 (stat) ± 0.16 (sys) GeV/c has been found to be in 
good agreement with the \J {p\) pair = 1-82 ± 0.09 GeV/c (see cq. 10). 

The analysis of two-particle correlation in Au + Au at ^snn= 200 GeV is also 
presented. Using the simple Glauber model we have demonstrated that the elliptic flow 
parameter v 2 cannot be generated solely by jet-quenching mechanism. The upper limit 
for v 2 from jet-quenching has been found to be less than 10%. It is evident that the 
large values of v 2 observed by all RHIC experiments have to be generated by other 
mechanisms. There is at least one model, inspired by "Color Glass Condensate" [14], 
which suggests the two-gluon correlations as a source of large azimuthal anisotropies. 
However; the elliptic anisotropy, in this model, does not respect the reaction plane. This 
seems to be in contradiction with observed correlation between the reaction plane and 
spatial distribution of particle emission after freeze-out observed by STAR collaboration 
[18]. An interesting alternative provides the "coalescence" model, where the partonic 
flow is "amplified" by means of parton recombinations [28] . This scenario is very exciting 
since it requires some kind of partonic collectivity in the initial state of the collision, 
which is anticipated if the quark-gluon plasma is created. 
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